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ABSTRACT 


An  attempt  was  made  to  model  the  transient  thermal 
response  of  jet  vanes  used  for  thrust  control.  A  simple 
computer  model  based  on  lumped  capacitance  methods  using 
boundary  layer  convection  and  stagnation  point  heating  as 
thermal  Inputs  appeared  to  adequately  predict  temperatures 
for  a  quarter-scale  model.  This  report  details  the  attempt 
to  enlarge  the  model  to  allow  comparison  between  thermal 
predictions  and  the  results  of  tests  on  a  full-scale 
prototype  jet  vane.  It  was  determined  that  the  model  could 
not  be  considered  a  thermal  representation  of  the  full-scale 
vane  assembly  and  several  modifioations  were  identified  in 
order  to  adapt  the  model  to  full-scale  applications. 
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I.  INTRODUCTION 


Because  of  the  need  for  maneuverability,  the  designs 
for  tactical  missiles  and  spacecraft  launch  vehicles  require 
the  application  of  active  control  systems.  Thrust  Vector 
Control  (TVC)  is  one  such  system  which  gives  the  capability 
of  trajectory  control  almost  independent  of  the  external 
forces  of  the  vehicle.  This  independence  is  important  when 
relative  air  flow  past  the  vehicle's  external  lifting 
surfaces  is  too  slow  to  generate  the  necessary  control. 
Applications  for  TVC  include  low  speed  flight  (during  launch 
or  while  hovering,  for  Instance)  and  some  cases  involving 
high  angle  of  attack.  Also,  in  cases  where  missiles  are 
launched  from  tubes  it  is  often  impossible  to  provide 
adequate  external  control  surfaces. 

One  method  of  TVC  is  the  insertion  of  a  jet  vane  into 
the  exhaust  of  a  rocket  nozzle  allowing  vehicle  control 
immediately  after  launch.  This  system,  while  allowing  for 
large  thrust  deflection  and  rapid  response,  also  results  in 
thermal  problems  and  some  thrust  loss.  Calculation  of  the 
heat  transfer  characteristics  is  made  difficult  due  to  the 
severe  thermal  environment.  So  a  reliable,  dynamic 
computational  model  would  be  beneficial  in  assisting  design 
efforts.  Such  a  model  of  the  thermal  process  in  the  jet 
vane  is  used  in  this  report.  This  model  was  developed  in 
support  of  a  larger  program  at  the  Naval  weapons  Center 
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(NWC) ,  China  Lake,  CA,  using  data  from  a  1/4  scale 
retractable  jet  vane  employed  in  the  stowable  Three-Axis 
Reaction  Steering  (STARS)  System,  Figure  1.1.  The  purpose 
of  this  report  is  to  verify  the  applicability  of  the  results 
gained  from  1/4  scale  testing  to  full  scale  prototypes. 


II.  BACKGROUND 


A.  THEORY 

The  traditional  approach  to  modeling  a  heat  transfer 
system  such  as  this  is  to  construct  a  comprehensive  model 
that  treats  the  flow  environment  of  the  vane  in  fine 
numerical  detail.  However,  a  complete  analysis  would  have 
to  account  for  multi-phase,  multi-component,  three- 
dimensional  and  time  dependent  effects  in  the  presence  of 
shocks,  boundary  layer  transition  and  turbulence,  separated 
flows,  surface  ablation,  chemical  reaction  and  solid-body 
and  gaseous  radiation  [Ref.  1].  Vast  simplification  can  be 
achieved  if  the  assumption  is  made  that  sufficient  accuracy 
can  be  gained  when  the  flow  and  vane  are  considered  to  be 
made  up  of  relatively  few  thermal  components  and  by 
realizing  that  the  net  effect  of  all  the  above  complications 
is  to  transfer  energy  to  and  from  the  vane  system.  The 
transient  thermal  responses  measured  at  accessible  locations 
on  the  vane  (protected  from  or  outside  the  exhaust  flow)  can 
then  be  used  to  estimate  temperatures  achieved  by  the  vane 
inside  the  flow. 

B.  THE  MODEL 

Work  on  the  jet  vane  thermal  model  was  begun  at  the 
Haval  Postgraduate  School  (NPS)  in  1986  by  Nunn  and  Kelleher 
[Ref.  1].  Development  of  the  model  was  continued  by  Nunn 


[Ref.  2]  and  Hatzenbuehler  [Ref.  3].  A  result  of  these 
studies  Is  a  model,  Figure  1.2,  of  the  jet  vane  undergoing 
testing  at  NHC.  The  important  features  of  this  model  are 
that  it  can  be  described  using  four  nodes  and  applying 
lumped  capacitance  procedures.  The  purpose  of  node  1,  the 
vane  tip,  is  to  account  for  stagnation  properties  of  thermal 
convection  near  the  vane  leading  edge.  The  fin  is  the 
remaining  portion  of  the  vane  exposed  to  the  rocket  exhaust- 
gas  temperatures,  and  is  modeled  by  node  2.  This  node  is 
subject  to  heat  transfer  by  both  turbulent  convection  and 
conduction.  Node  3  is  the  vane  shaft  with  heat  transfer 
assumed  to  take  place  by  conduction  only  and  node  4 
represents  the  mount  connecting  the  shaft  to  the  rocket 
frame,  also  subject  to  conduction  only.  Details  of  the 
model  are  contained  in  Reference  2,  pages  9-11,  and 
Reference  3.  The  governing  equations  for  heat  transfer  in 
this  model,  become  [Ref,  2,  p. 31,  Ref  3]: 

T,  -  -  (T,/C,)(l/R„  f  l/R,j)  +  Tj(l/C,R,j)  (1) 

+  T„(1/C,R„) 

Tj  “  (VCjXl/R^j)  -  (Tj/Cj)(l/R,j  +  1/R,j  +  i/R^)  ^2) 

+  (VC2)(1/R25)  *  (T^/CjXl/R^j) 

Tj  (Tj/Cjyd/Ra)  “  (VCs)(VRa  ♦  i/Rj*)  (3) 

+  (VCsXVR,*) 

T*  “  -  (VCJd/RsiWRui)* 
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Figure  1.2.  Four  Kode  Model  Configuration 
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Here,  the  subscript  R  refers  to  recovery  temperatures  which 
are  used  because,  due  to  the  high  speed  in  the  compressible 
boundary  layer,  the  temperature  driving  the  heat  transfer 
will  exceed  the  local  static  temperature  of  the  exhaust.  It 
>  should  be  noted  that  each  of  the  temperature  coefficients  on 

the  right-hand  sides  of  Eqs.  (l)-(4)  have  the  dimensions  of 
inverse  time.  In  fact,  the  RC  products  are  representations 
of  the  time  constants  describing  the  energy  transport 
processes  occurring  at  and  around  the  nodes.  Among  the 
several  resistances  and  capacitances  appearing  in  the  model, 
those  believed  to  be  roost  uncertain  in  their  values  are  the 
input  resistances  due  to  convection,  and  the  time  constants 
associated  with  conduction  through  the  mount. 

^  System  identification  procedures  were  used  to  determine 

what  these  values  should  be  in  order  to  provide  agreement 
'  with  temperature  histories  recorded  for  the  shaft  and  mount 

during  NWC  test  firings  on  the  1/4  scale  vanes.  The  Initial 
values  used  for  resistance  and  capacitance  are  listed  in 
Appendix  A. 

C.  SYSTEM  BUILD 

It  was  decided  to  construct  the  simulation  model  using  a 
personal  computer  (PC)  version  of  System  Build  and  System 
Identification  developed  by  Integrated  Systems,  Inc. 

»  [Ref.  4,  5].  In  order  to  simplify  the  c<.)de,  the  following 

parameters  were  defined  for  use  in  Eqs.  (l}-(4): 
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aai  “ 

VC2R,2 

aaa  =  aj,  +  a23  +  bjj 

^23 

—  I/C2R23 

a32  “ 

I/C3R23 

333  =  (^32  +  334) 

^34 

=  I/C3R34 

11 

VC.Rj, 

^44  = 

^40 

*  I/C4R4G 

bii  = 

VCiR,, 

baa  ~  ^/^2^F2 

These  are  referred  to  as  the  characteristic  rates  of  the 
model,  and  the  resulting  equations  for  input  to  system  build 
become : 


T,  = 


“  ^22^2  ■*'  ^23^3  ^22'^R2 


-  a33T3  +  aj^T^ 


a43'^3  "  a^^T^ 


(5) 

(6) 

(7) 

(8) 


The  System  Build  thermal  model  consists  of  nine  Super 
Blocks : 

NODI IN,  N002XN,  N0D3IN,  N0D4IN, 

NODEl,  N0DE2,  NOOE3,  N0DE4, 

VANE, 

as  shown  in  Appendix  B.  The  input  is  a  ramp-up,  plateau, 
ramp-down  to  provide  for  the  transient  stagnation  tempera¬ 
ture  profile  generated  by  the  rocket  motor  firing,  [Ref.  2, 
p.23].  The  first  four  blocks  compute  the  parameters 
necessary  for  the  equations.  The  ramp-up,  plateau,  ramp- 
down  input  to  NODllN  and  N0D2IN  has  a  maximum  value  of 
unity.  In  NODlIN,  for  example,  this  input  is  multiplied  by 
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and  as  gains  to  form  the  time  varying  element  in 
eg. (5).  Also,  in  NODIIN,  a^2  generated  by  a  step 
function  and  added  to  b,,  to  form  a„.  So  the  outputs  of 
NODIIN  are  a,,,  and  a^j* 

In  NODEl  the  outputs  of  NODIIN  are  combined  with  the 
external  input,  Tj,  to  form  eg. (5).  The  integrator  then 
converts  to  which  is  the  NODEl  output.  N0DE2,  NODE3 
and  N0DE4  are  similar  in  form  and  function  and  generate  Tj, 
Tj  and  respectively.  The  super  block,  VANE,  inter¬ 
connects  these  last  four  blocks  and  provides  the 
simultaneous  solution  for  the  four  nodal  temperatures.  The 
step  amplitude  and  gain  parameters  in  any  of  the  blocks  may 
be  varied  to  facilitate  system  identification.  When  a 
parameter  is  allowed  to  vary,  System  ID  will  calculate  the 
value  which  allows  the  closest  fit  of  the  model  predictions 
to  actual  test  data  curves. 

This  simulation  was  run  initially  with  only  a  3-node 
model  using  the  mount  as  ground  for  reference  [Ref.  2]  and 
setting  the  mount  thermal  resistance  as  an  adjustable 
parameter.  Because  of  a  fault  in  connecting  the  System 
Build  super  blocks,  erroneous  results  were  obtained  for  the 
values  of  b22  and  a^)  [Ref.  2,  p.  453.  After  removing  the 
fault,  the  following  corrected  values  were  obtained: 
bjj  •  0.19S6  s'^  and  a^j  -  -0.4627  s*^.  These  values  were 
obtained  through  System  Identification  using  the  NWC 
transfer  function  for  Tj  [Ref  2,  P.  24]  for  comparison  with 
the  simulation  results. 
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III.  APPROACH 


A.  STARTING  POINT 

The  coefficient  values  obtained  by  Hatzenbuehler 
[Ref.  3]  were  also  based  on  faulty  System  Build  model,  and 
these  numbers  were  re-calculated  before  attempting  to  scale 
up  the  model.  Using  event  2  data  from  Reference  3 
[Appendix  C] ,  the  model  was  again  set  up  for  system 
identification  with  the  MAXLIKE  function  of  MATRIXx  [Ref  5]. 
Four  parameters  were  allowed  to  vary  simultaneously  and  were 
compared  to  test  data  temperatures  recorded  from  the  shaft 
(Node  3),  and  the  mount  (Node  4).  Initial  values  for  all 
parameters  are  contained  in  Appendix  D.  The  parameters 
which  a  were  allowed  to  vary  were: 

b;.  -  allows  estimation  of  the  value  of  R^jr  the 

resistance  to  heat  transfer  to  the  fin  surface, 

a34  -  allows  estimation  of  R34,  the  resistance  to 
conduction  through  the  shaft, 

a^;  -  given  a^,  allows  estimation  of  the  capacitance  of 
the  mou^,  C^,  and 

a^Q  -  allows  an  estimate  of  the  conductive  resistance  to 
hea^  cransfer  through  the  mount  to  the  environment, 

Three  (3)  iterations  of  MAXLIKE  returned  the  values  of: 
b22  -  0.2718  S*’,  a34  »  0.299  S  a^j  -  0.1322  s’ 

and 

a^g  *  O.lOlO  b'\ 
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These  correspond  to  values  of: 

Rf2  =  1.67  K/W,  =  3.33K/W, 

=2.27  J/K  and  R4Q  =  4.33  K/W. 

Figure  3.1  shows  the  excellent  fit  these  values  give  to  the 
test  data.  It  should  be  noted  that  because  of  these  new 
maximum-likelihood  parameter  values,  the  simulation  results 
shown  here  are  slightly  different  from  those  reported  in 
Reference  3.  For  instance,  the  maximum  vane  tip  temperature 
predicted  for  event  2  is  1438  K  (1144  K  above  ambient) , 
some  175  K  higher  than  that  reported  in  Reference  3. 

B.  ENLARGING  TO  FULL  SCALE 

In  the  time  since  this  study  began,  NWC  China  Lake  has 
begun  testing  of  full-scale  prototype  jet  vanes.  With  the 
larger  vane  it  became  possible  to  insert  thermocouples 
inside  the  fin  in  locations  corresponding  to  some  of  the 
nodes  in  the  thermal  model.  The  only  node  still 
Inaccessible  for  actual  measurement  was  the  vane  tip,  where 
the  highest  temperatures  were  expected  to  occur.  Therefore, 
an  attempt  was  made  to  "scale  up"  the  computer  thermal  model 
to  compare  with  the  results  of  the  prototype  testing. 

Three  important  differences  between  the  sub-scale  and 
full  scale  testing  configurations  must  be  noted.  In  the 
1/4-scale  testing,  four  fixed  vanes  were  arranged 
symmetrically  in  the  jet  nozzle,  only  three  vanes  were  used 
in  the  full-scale  firing  and  one  of  these  was  deflected 
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Figure  3.1.  Temperature  Histories,  Pour  Node  Model 


intermittently  during  the  test.  Though  these  vanes  were 
also  arranged  symmetrically,  it  is  likely  that  the  heat- 
^  transfer  rates  were  significantly  affected  by  the  different 

shock-wave  patterns  generated  through  the  vane  arrangement. 
Additionally,  although  the  full-scale  vanes  were  assumed  to 
have  been  strictly  proportional  to  the  sub-scale  model,  the 
mounting  apparatus  was  not.  In  fact,  no  consistent  mount 
arrangement  was  used  during  the  1/4-scale  testing  [Ref  3]. 
Since  consistent  data  for  the  resistance  to  heat  transfer  in 
the  mount,  and  from  there  to  ground,  was  unavailable,  it  was 
decided  to  return  to  the  3-node  thermal  model  [Ref  2]  which 
connects  the  shaft  node  directly  to  ground.  This  model  is 
summarized  in  Appendix  E.  Without  geometric  similarity 
I  between  the  two  mounting  configurations  there  was  little 

reason  to  believe  that  the  time  constants  relating  to  the 
^  mount  could  be  scaled  directly.  Finally,  in  video 

recordings  [Ref.  6]  of  the  full-scale  tests,  both  the  vane 
mounts  and  the  jet  nozzle  wall  exhibited  Indications  of  more 
rapid  heating  than  in  the  subscale  tests.  The  prototype 
mounting  apparatus  included  fin-like  extensions  parallel  to 
the  flow  and  subject  to  convective  heat  transfer.  Also,  the 
nozzle  wall  was  proportionally  much  thinner  in  the  full- 
scale  firings.  This  would  have  allowed  the  nozzle  to  reach 
higher  temperatures.  As  will  be  discussed,  this  work  has 
shown  that  these  effects  may  have  significantly  altered  the 
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radiation  environment  of  the  jet  vanes  resulting  in  faster 
heat-up  transients  and  a  longer  cool-down  period  following 
motor  burnout. 

For  the  characteristic  coefficients  affecting  only  the 

vane  and  shaft  nodes,  it  was  felt  that  direct  scaling  was 

appropriate.  The  characteristic  rate,  a^j/  in  the  1/4  scale 

model,  for  example,  is  decreased  to  its  full-scale  value  by 

the  following  method: 

VC,R,j  p  -  density 

C  =  p  CpV  Cp  -  specific  heat 

R  »  L/k  A  V  -  volume 

L  -  Length 

k  -  thermal 

conductivity 

A  -  area 

Let  the  subscript  f  denote  full  scale  values  and  s  the 
values  corresponding  to  reduced  scale.  Then  knowing  that 
the  thermal  diffuslvity,  a«k/pCp,  remains  constant  during 
scaling: 

ai2,  »  o.V(V,  L.) 

-  a  A^  (scale)V[Vf  (scale)^  L|  (scale)] 

or 

(scale) ^ 

In  this  case  (scale)  is  1/4  and  a,2t'°^*‘^  ^  so 
“  (4.7)  (1/4)*  »  0.294  S*’. 
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The  remainder  of  the  characteristic  rates  governing 
conduction  in  the  model  were  similarly  scaled. 

The  characteristic  rates  governing  convection  were 
scaled  according  to  the  procedure  outlined  in  Reference  2 
(p.20).  Stagnation  point  thermal  resistance  decreases  as 
(0.25)*’’®.  However,  upon  review  of  the  equations,  it  was 
determined  that  boundary  layer  thermal  resistance  should  be 
proportional  to  (0.25)'^*®.  For  instance,  the  stagnation 
point  characteristic  rate  is  given  by; 
b,,3  =0.45  s’’ 

«  1/[C^  (0.25)®R„  (0.25)*'-®] 

=  1/C^R^^(0.25)'''® 

-  b„/(0.25)’*® 

and, 

b,„  -  (0.25)’-®(b^„)  »  0,0563  s"’. 

The  entire  set  of  coefficient  values  for  the  full-scale 
thermal  model  is  given  in  Appendix  F. 

C.  SIMULATION 

Additional  Information  required  before  the  computer 
model  can  be  executed  includes  the  recovery  temperatures  and 
the  input  function  modeling  the  thrust.  From  evaluation  of 
the  propellant  properties,  NNC,  China  Lake  calculated  the 
stagnation  temperature  as  being  2971  X.  Ambient  temperature 
was  not  reported  but,  based  on  the  initial  readings  of  the 
21  thermocouples  used,  and  accounting  for  their  differences 
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according  to  their  location,  was  chosen  as  301  K. 
was  then  calculated  using  the  procedures  of  Reference  2 
(pp.  17-21) .  Referring  the  recovery  temperatures  to  ambient 
yields  =  2670  K  and  1^2  =  2567  K.  These  are  the  values 
used  in  the  model. 

The  thermal  input  schedule  was  modified  to  synchronize 
with  the  thrust  trace  obtained  from  the  full  scale  firing  as 
shown  in  Figure  3.2  The  resulting  normalized  thrust 
function,  for  input  to  super-blocks  NODIIN  and  N0D2IN  of  the 
model  is  given  in  Table  I. 

As  can  be  seen  from  Figure  3.3,  simulation  with  the 
model  under  the  assumption  of  full  thermal  scaling  failed  to 
predict  the  temperature  of  node  2,  much  less  to  give  a 
reasonable  estimate  for  the  temperature  of  node  1,  the  node 
of  interest*  This  leads  to  the  conclusion  that  the 
differences  between  the  two  (large  and  small)  devices  are 
not  adequately  described  by  a  single  scale  factor.  However, 
investigation  proceeded  to  determine  If  the  model  could,  in 
fact,  still  serve  the  purpose  for  which  it  was  constructed: 
to  predict  the  maximum  temperatures  at  the  vane  tip  and 
other  locations  of  special  interest. 


Figure  3.2.  IMC  Thruat  Trace,  Prototype  Testing 


TABLE  I. 

Input  Thrust  Function 
TIME(S)  THRUST  (%) 


0 

0 

0.07 

81 

0.14 

78 

1.67 

100 

5.81 

100 

6.0 

94 

6.35 

94 

6.70 

0 

120.0 

0 

18 


19 


Figur*  3.5.  Tenperature  Histories  Based  on  Conventional  Scaling  Procedures 


IV.  REVISED  PROCEDURE 


A.  EXAMINATION 

In  order  to  gain  further  insights  from  the  model  it  was 
first  necessary  to  realize  where  the  probable  inconsis¬ 
tencies  occurred  in  the  assumption  that  it  could  be  a  scaled 
representation.  One  source  of  possible  error  is  related  to 
the  fact  that  in  modeling  the  subscale  process  [Ref  2,  3] 
the  net  theirmal  input  resistances  were  increased  to  values 
well  above  those  based  solely  upon  thermal  convection.  Such 
increases  were  found  to  be  necessary  to  obtain  good  agree¬ 
ment  with  the  1/4 -scale  data,  but  the  present  model  does  not 
provide  a  rationale  for  adjusting  these  factors  on  the  basis 
of  scale.  Such  a  rationale  would  have  to  come  from  a 
thermal  model  that  includes  the  effects  of  radiation  and 
ablation  on  the  net  thermal  input  resistance.  Further,  as 
previously  mentioned,  the  radiation  effects  in  the  two 
testing  configurations  appeared  to  be  drastically  different. 

In  addition  to  the  input  resistance,  the  thermal  heat 
sink  effect  of  the  mount  (output  resistance)  was  also  not 
amenable  to  scaling.  Since  variations  of  these  factors  can 
have  significant  impact  on  the  agreement  of  the  model  with 
known  data,  it  was  decided  to  apply  System  Identification 
with  the  known  shaft  temperature  of  the  full  scale  prototype 
to  see  if  matching  this  temperature  would  successfully 
predict  the  temperatures  at  nodes  elsewhere  on  the  test 


vanes.  <^nce  again  parameters  ajj  and  b22  were  allowed  to 
vary.  Full  scale  values  [Appendix  F]  of  the  other 
parameters  were  used. 

B.  RESULTS 

Figure  4.1  shows  the  lack  of  success  with  this  approach. 
While  allowing  the  modeled  temperatures  at  node  3  to 
approach  the  data  for  the  shaft,  the  remaining  temperatures 
increased  disproportionately,  to  the  point  where  the  node  2 
•*:emperature  reached  its  maximum  value,  What  became 

apparent  was  that  the  temperature  profiles  of  the  data  for 
the  full  scale  testing  were  quite  different  from  those  for 
the  subscale  test  firings.  The  rapid  increase  in  the  data 
temperatures  in  the  first  ten  seconds  more  closely  followed 
the  expected  rise  for  those  nodes  receiving  direct  thermal 
input  from  the  surroundings.  Thus,  it  became  apparent  that 
in  the  full-scale  tests,  the  shaft  was  receiving  significant 
thermal  input  in  addition  to  that  due  to  conduction  through 
the  vane. 

C.  CONCLUSION 

Because  of  the  change  in  testing  conditions  the  1/4- 
scale  model  cannot  be  used  directly  to  predict  the  required 
temperatures.  Some  source  of  direct  heat  input  to  the  shaft 
(e.g.,  via  radiation)  seems  to  be  indicated.  Without  this 
direct  heat  input,  any  attempt  to  match  the  data  temperature 
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Figure  4,1.  Temperature  Histories  A.f ter  SYSTEM  ID 


for  the  shaft  to  the  model's  predictions  causes  vast 
distortion  of  predicted  temperatures  elsewhere  on  the  vane. 
The  two  test  configurations — "large"  and  "small" — were  not 
scale  models  of  each  other  from  a  heat  transfer  point  of 
view. 


APPENDIX  A 
[Reference  2] 

Resistances  and  capacitances  for  1/4 -seals  model 

Rf,  =  20.3  K/W 
Rp2  =  2.3  K/W 
*R23  =  2.1  K/W 
C,  =  0.11  J/K 
C2  =  2.2  J/K 
C3  =  1.0  J/K 

*R23  is  revised  to  agree  with  the  corrected  3 -node  model 
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APPENDIX  B 


) 

) 


System  Build  Blocks  For  4  Node  Model 


) 
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APPENDIX  C 
[Reference  3  APP.C] 

NWC  Test  Firing  Parameters 


) 


) 


Propellant 
Molecular  Weight 
K= 

C*= 

Burn  rate  Coefficient  = 
Burn  rate  exponent  = 
Ambient  Conditions  = 


0%  AL  HTPB 
26.1 
1.210 

4930  ft/sec 
0.934 

.05703  in/sec 
.28063 

70  F,  13.7  PSIA 


Propellant  Mass  = 
Propellant  Density  ® 
Initial  Throat  Area  = 
Post-fire  Throat  Area  =» 
Max  Pressure  = 

Max  Thrust  ® 

Total  Impulse  a 
ISP  = 

Exit  Dia  » 

Max  Exit  Pressure  » 


C.17  Ibm 
.06185  Ibm/in* 
.14750  in* 
.15135  in* 

2266.8  PSI 
514  Ibf 
1535  lb-sec 

248.8  sec 
1.8125'* 

12.3  PSI 


t 
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APPENDIX  D 

Initial  Values  For  1/4  Scale  Model 
NODIIN 

A12  =  4.7  s*’ 

Bll  =  0.45  &'' 

TRl  =  2360  K 

N0D2IN 

A21  =  0.23  s'^ 

A23  =  0.18  s‘^ 

*B22  =  0.1956  s''' 

TR2  =  2260  K 

NOD3IN 
A32  =  0.38 
*A34  »  0.17 

N0D4IN 
*A43  »  0.17 
*A4G  -  0,14  ; 

♦paraneter  allowed  to  vary  for  SYSTEM  ID 
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APPENDIX  E 
[Reference  2] 

Three  Node  Heat  Transfer  Model 


6 


Figure  S. 1 .  Basic  Three  Node  Model 
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Governing  equations 

T,  =  -  (VCi)(l/Rf,  +  I/R12)  +  (T2/C2)  (1/R,2)  +  WqRfi 

T2  =  (T,/C2)(1/Ri2)  -  (T2/C2)  (1/Rf2  +  I/R12  +  VR23)  +  VC2R23 
+  (WC2)  (1/Rf2) 

Tj  =  (T2/C3)  (I/R23)  -  (T3/C3)(1/R23  +  1/R3e) 

Characteristic  Rates: 
a,i  =  a, 2  + 

^12  “  VC,R,2 
^21  ~  VC2R12 

®22  “  ^21  ®23  ^22 

®23  “  VC2R23 

^32  “  ^23 

^33  “  ^32  +  b33 

b^^  =*  1/C^Rf^ 

^22  “  ^/^2^F2 

^33  *  VC3R36 

Equations  for  System  Build 
Ti  =>  -a,,T,  +  a,2T2  +  b^, 

Tj  «  aj^T,  -  a22T2  +  ajjTj  +  b22 
T3  “  ^32^2  -  a33T3 


V 
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APPENDIX  F 

Full  Scale  Parameters 


Characteristic  Rates: 

a, 2  =  0.294 
aji  =  0.0144 
agj  =  0.01125 
ajj  =  0.0238 
ajj  =  0.030 
=  0.0563 
b22  =  0.037 

Full  Scale  Temperatures: 


Tri  =  2670  K 
Tg2  =  2567  K 
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